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% Justification de choix de sujet : 


J'ai porté la réflexion sur ce sujet car il est fréquemment rencontré en génie civil et 
d'une importance primordiale. L'usage des premiers murs de soutènements remonte à 
l'antiquité principalement dans les monuments comme les amphithéâtres et joue le rôle de 
gradin des talus sur les collines où se situait l'ouvrage, ensuite dans le secteur agricole pour éviter 
l'érosion de sol avec des terrains pierreux en pente. 


Dans l’époque moderne, et sous l'effet de l'accroissement démographique et l'exode urbain la 
méthode s’est ensuite étendue pour protéger les constructions d’un éventuel glissement ou un 
effondrement de terrain ou préserver les routes et chemins des éboulements. De ce fait l'étude 
de cet ouvrage doit viser la sécurité, moins dépenser de la matière et économiser le cout de la 


réalisation. 





image 21: TEE encre d’un amphithéâtre 


de l'empire Romaine Image 2 : Culture en terrasse 





Image 3 : Mur de soutènement sur un terrain accidenté 
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% Introduction : 


- Mur cantilever est un mur en béton armé reconnu par sa forme de T renversé ou L 
caractérisé par un voile encastré dans la semelle de fondation et une base élargie encastrée 
dans la partie supérieure du sol afin de faire participer le remblai à soutenir dans sa stabilité, 
et ci-dessous des images lors de sa construction pour un projet ferroviaire et routier dans 
zone montagneuse. 


. :semelle 


eu 





- La construction doit satisfaire les conditions géotechniques et structurelles pour assurer sa 
résistance, sa durée de vie et son utilisation en service, car il fait partie du réseau routier et 
permet la sécurité des usagers sur la route et donc éviter de nombreux accidents. 


Toutefois Le volet économique est aussi important et fait partie de la conception de ce 
fait l'ingénieur doit réaliser une conception appropriée qui doit satisfaire à toutes les 
exigences et fait l’objet de nombreuses études, où l'optimisation intervient en essayant de 
traiter des problèmes par exemple : diminuer le poids, la taille et également le coût de 
production et les résoudre avec des méthodes d’optimisations. 


> Éclaircissement de la Problématique : 


La conception d’un mur de soutènement était basée au paravent sur des méthodes 

classiques tel que le principe de l'essai et erreur, on propose des dimensions pour les 

parties formant l'ouvrage en fonction de la hauteur du remblai et en satisfaisant aux 

exigences géotechniques vis-à-vis du: 
Le Rupture 


Le Renversement 
Le Poinçonnement. 
+ 


+? Glissement a) Poinçonnement b) Renversement c) Glissement 
Si cette conception n’est pas vérifiée A 


un redimensionnement de l'ouvrage 
sera nécessaire 





d) Rupture de la structure e) Instabilité générale du site 


Donc pour avoir une conception optimale vis-à-vis de coût et sécurité n’est pas sûr 
certain et refaire le calcul constitue une perte de gain et de temps, pour cette raison, 
intervient l'utilité des méthodes d'intelligence artificielle. 


Diverses méthodes de modélisation mathématique peuvent faire cette conception à 
coût minimal à savoir : l'algorithme génétique, optimisation par essaims particulaires 
ou algorithme de colonie de fourmis..., ils s’utilisent pour résoudre des problèmes non 
linéaires. 


Plus particulièrement on choisit d'employer l'algorithme génétique qui fait partie 
des méthodes métaheuristiques permettant d'assurer une conception optimale 
économique et sécurisante pour l'ouvrage, qui s'inspire de la théorie de l’évolution de 
Darwin en utilisant une analogie entre un problème génétique et un problème 
d'optimisation. 


Dans ce qui suit on va traiter le principe de fonctionnement de cette méthode et 
son application sur un mur de soutènement cantilever. 


> Explication du Principe de fonctionnement : 
Définition de la méthode métaheuristique : 


- Le mot métaheuristique est dérivé de la composition de deux mots grecs : "heuristique” 
qui désigne trouver et le qualificatif "méta" qui signifie adaptable à un nombre de problèmes 
sans changements majeurs dans l'algorithme. 

- C'est un même schéma de recherche appliqué plusieurs fois au cours de l'optimisation, et 
directes, en utilisant une population de points comme méthode de recherche et souvent 
inspirée par des systèmes naturels, qu'ils soient pris en physique ,ou en biologie de 
l’évolution tel que les algorithmes génétiques. 


Définition de l'algorithme génétique : 


C'est un algorithme d'exploration et inspiré du concept de sélection naturelle de Darwin et 
a été proposé par Holland en 1975 pour décrie les systèmes adaptatifs, le concept de cette 
méthode est de générer une descendance, au début la population est soumise aux 
contraintes du milieu naturel et seul les éléments les mieux adaptés survivent, ensuite de 
nouvelles générations apparaissent et la sélection naturelle permet de trier les parties des 
meilleurs éléments de la génération précédente, en utilisant des règles de transition 
probabilistes afin de guider l'exploration 


le vocabulaire technique employé est lié à la théorie l’évolution de de la génétique : La 
solutions codée —chromosome, variables de décision— (gène) tandis que les points de 


l'espace de recherche sont alors les individus d’une population et la fonction à optimiser 
correspond à leur adaptation (fitness). 


Espace de recherche 


Population 


Individu 


Schématisation des niveaux d'organisation d’un algorithme génétique. 





Afin d’expliciter le fonctionnement de cette méthode d'optimisation, sur la figure une 
représentation des différentes étapes d’un algorithme génétique simple : 





Initialisation de la 
population 






Reproduction 


(croisement et mutation) Évaluation des 


individus 


oo Non 
Phase sélection o r 
Critère d’arrêt 





Oui 


Détermination des 


solutions 





Explication de schéma illustrée : 


1)- Initialisation: l'algorithme démarre avec une population initiale d'individus, le 
choix des individus conditionne fortement la rapidité de processus. Si la position de 
l optimum dans l’espace de recherche est totalement inconnue, il est intéressant que 
la population soit répartie sur tout l’espace de recherche. 


2)- Évaluation des individus : Consiste à évaluer chaque solution contenue dans la 
population, la mesure de performance des solutions s'appuie sur la valeur de la 
fonction objective, cette étape permet de classer les solutions, afin de déterminer les 
solutions qui seront sélectionnées pour construire une nouvelle population de 
solutions. 


Fonction objectif: c’est la fonction de coût ou “fitness” dans la communauté 
d'’algorithmes évolutionnaires qui définit le but à atteindre, souvent dans les 
problèmes d'optimisation on chercher à minimiser la fonction objective f (x), tout en 
respectant les contraintes qui peuvent être linéaires ou non linéaires et illustré ci- 


dessous : 
gi(x) <0 iell:p] Avec: dix) : contraintes d'inégalité 
h;j(x)=0 jeliml] h;(x) : contraintes d'égalité 
Lk <X, SU, kEfl:n] Lk, U, : contraintes aux limites. 


GS cc 
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3)- Stratégie de sélection : Pour objectif d'identifier les individus qui doivent se 
reproduire sur la base de leur fonction d'adaptation, les individus les mieux adaptés 
sont sélectionnés alors que les moins bien adaptés sont écartés, Il existe plusieurs 


techniques de sélection tel que : sélection par roulette, sélection par tournoi et 
sélection par rang. 


i) Sélection par roulette : la solution obtenue est copiée dans la nouvelle population 
proportionnellement à sa fitness, en effectuant des tirages avec remise pour les 


individus existant dans la population, la probabilité de sa sélection dans la nouvelle 
f(ch;) Avec : 


Zi=1 (Chi) ch; : individu 
f(ch;): fonction d'adaptation 


population de taille N est: P(ch;) = 





séléction pa rouelette. 


ii) Sélection par tournoi : Le concept est de tirer deux individus aléatoirement dans 
la population et on reproduit le meilleur des deux dans la nouvelle population. On 
répète la procédure jusqu’à ce que la nouvelle population soit complète. 


iii) Sélection par rang : C’est une procédure similaire à celle par roulette et consiste 
à classer la population suivant la fonction d'adaptation, chaque individu se voit 
accorder un rang. Plus l'individu est meilleur, plus son rang est élevé. L'ensemble des 
individus est représenté sur un segment de droite dont les valeurs sont comprises 
entre O et 1. 


4)-Reproduction : Une fois que le choix des individus est effectué, le rôle de la phase 
de reproduction est l'application d'opérateurs de variation tels que la mutation et le 
croisement. 


i) Opérateur de mutation : Le rôle de cet opérateur est de changer aléatoirement la 
valeur d’un individu. Dans le cas du codage binaire, chaque bit a; € {0, 1} est remplacé 
par son complémentaire a; = 1-a; . Dans l’image ci-dessous, une mutation a eu lieu 
sur le cinquième gène du chromosome et elle a transformé de 1 en 0. 





héma de mutatior 


Cet opérateur est lié par deux hypothèses : 

- le gène peut prendre des valeurs dans un intervalle xn, ao, 
- [a probabilité de mutation pm est fixée. 
Après on procède au tirage au sort pour chaque gène, et Il est de deux façon : 
- [la mutation uniforme : consiste à remplacer le gène x; sélectionné par une valeur 
quelconque x’, tirée aléatoirement dans l'intervalle [a Lg À 
- [a mutation non uniforme : Le gène x, subit des modifications durant les premières 
générations, puis graduellement décroissantes au fur et à mesure que l’on progresse 
dans le processus d'optimisation. 
Donc pour une génération t, un tirage au sort est fait sur une valeur binaire qui 
décidera si le changement doit être positif ou négatif, la nouvelle valeur x’, du gène 
Xg + A(t, XP — xg) si rand = 1 

k — A(t, Xg — 0) sinon 
Avec : A(t, y) est une fonction qui définit l'écart entre la nouvelle valeur et la valeur 
initiale à la génération t 
- rand : un nombre aléatoire qui prend les valeurs O ou 1 
ii) Opérateur de croisement : cherche à créer des nouveaux individus par l'échange 
d'information entre les chromosomes avec le biais de leur combinaison, ensuite la 
population est divisée en deux sous populations en conservant la taille, après le couple 
est formé par un membre provenant de chaque sous population et participe au 
croisement avec une probabilité supérieure à 50%. 
Si le croisement a eu lieu entre les chromosomes parents (ch, et ch), on tire 
aléatoirement une position de chacun des parents. 


x, est donnée par: x'g 


eı| | [T|] 


N # 


x > a — ET" 
© CROISEMENT?) 
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FIGURE 2.4 — Croisement standard en un seul point 





D'où l'échange se fait entre les deux sous chaînes terminales de chacun des 
chromosomes, ce qui produit deux enfants ch’, et ch', comme illustrée sur l’image 
ci-dessous. 

- Une remarque très importante : le croisement en un seul point prévue à la quatrième 
position a eu lieu entre les chromosomes ch, et ch, . 

Le processus de combinaison peut s'étendre aussi en multipoints indiqué sur l’image 
ci-dessous : 


—— 


EE Mie ~ 
( CROISEMENT) 


A 


FIGURE 2.5 — Croisement standard en deux points 





Ce croisement est très efficace et peut s'étendre à n'importe quel type de chaînes 
(réelles, binaire, ..) 

Apparemment dans le cas des chaînes réelles le croisement de type barycentre est 
préférable et consiste à sélectionner deux gènes ch1(i) et ch2(i) dans chacun des 
parents à la même position i pour créer deux nouveaux gènes ch’ 1 (i) et ch’ 2 (i) par 
combinaison linéaire 


cha(i) = (1 — a) x chi (i) +a x cho(i) 


{ chili) = a x chili) + (1 — a) x chali) 


- L'étape suivante consiste à déterminer les principaux paramètres pour évaluer la 
conception optimale de cet ouvrage en satisfaisant les exigences de 
dimensionnement. 


-Après on cherche une combinaison de ces paramètres (coût, spécifications et 
fonction) afin de suivre une procédure de maximisation ou de minimisation. 


- Afin de résoudre un problème d'optimisation avec contraintes nécessite souvent la 
programmation de la méthode de calcul avec un logiciel de programmation avancé 
pour optimiser le temps de recherche, pour cela on va créer un fichier contenant les 
données de l'algorithme puis les implémenter en MATLAB. 
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% Application de l'algorithme génétique pour l'optimisation d’un mur 


cantilever : 





Fiqure 1 
La figure 1 ci-dessus montre un mur de soutènement modélisé par 5 variables de 


conception et soumis à sept forces : 
W.: le poids combiné de toutes les sections du béton 


armé mur. 
W;: le poids du remblai agissant sur le talon du mur. 
W: le poids du sol sur le pied du mur. 


- X. :Largeur de la semelle 
- X- :Largeur de la pointe 


- X> :Épaisseur de voile Q : une surcharge 
- À, ‘Epaisseur en haut de voile P, : Ja force résultante de la pression latérale active. 
- Xs ‘Épaisseur de la semelle P, : la force résultante de la pression passive de la terre 
P, : la force résultante de la contrainte portante du sol 
de base. 


Formulation mathématique du problème : 


- Pour atteindre la conception optimale de l’ouvrage on définit la fonction objective f (x) 
qui est le volume de l'ouvrage, et selon la figure ci-dessus le volume du mur peut être 
donné par l'équation : 


(X3 + À5) 
FO = X x X, + xH En 
H : la hauteur de mur 
On remarque que notre problème est mono-objectif, et le but c'est de minimiser la 
fonction objective f (x) ou simplement le volume de mur tout en respectant les contraintes 


de stabilité de la géométrie et les exigences de dimensionnement. 
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c- Contraintes 


1.Contraintes de stabilité (Théorie de Rankine ) Voir annexe: 


2.Contraintes géométriques : 
- Les contraintes géométriques dirigent l'algorithme à chercher des solutions 


optimales dans un certain intervalle fixé par des limites inférieures et supérieures de 
chaque partie de l'ouvrage, et pour faire des conceptions convenables sont 
constituées des contraintes aux limites, et ci-dessous une formulation tirée de la figure 
1 sur les contraintes : 

- Xi =X- +X +X6 

- X3> X; 

- Xę> À) 

- X6> X3 
- Les équations ci-dessus représentent des contraintes linéaires, d'égalité ou 
d'inégalité. 
- La méthodologie algorithmique adoptée dans cette partie est présentée par 5 
étapes : 


Étape 0 : 











Donnes initiales : 
-Données du sol remblai et 
sol de fondation 
-Données du mur (poids 
volumique) 
-Le nombre de variables 
(Dimensions du mur) 
-Domaine de variation des 
variables 

-Nombre de population 

-Critère d'arrêt 


Étape 1 : 










Initialisation d’une 
population initiale 
aléatoire 













Étape 2 : 



















Vérification de la 


stabilité en calculant 
les facteurs de sécurité 





Étape 5 : 


Étape 4 : 


Détermination de la 
solution qui donne le 
Volume minimum (Avec 
logiciel MATLAB) 


Répétition du 
processus 


Étape 3 : 






algorithmique jusqu’à 
arriver à l'étape 
d'arrêt 





Évaluation de 
volume de l'ouvrage 
pour chaque 
solution trouvé 






Ci-dessous un exemple de validation appliqué dans un projet de fin d’études en 
s'aidant d’un logiciel de programmation MATLAB : 
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d-Application de Processus algorithmique sur L'ouvrage Cantilever : 

Objectif : montrer l'efficacité et la performance de l'algorithme génétique dans l'optimisation de 
volume de mur. 

- C'est un mur d’une hauteur de 5,65m et un angle de remblais a°=0° et d’une surcharge 

Q= 10,2kN/m2. 

- Ci-dessous les valeurs des paramètres utilisés : 


Les Paramètre de Sol deremblai Sol de Fondation fco(MPa) 25 
sol 


A] RE U(coefficient de 0,3 M 
poisson) 
C T N/m?) = 
Tableau 2 : Les Caractéristiques de 


Tableau 1 : Les Paramètres de sol materiaux 





Avec : 
fe : résistance à la compression du béton.  E;;: module de Young 


C : Cohésion du terrain 
y :Poids spécifique total de sol utilisé 


-À l’aide de Logiciel Matlab l’évolution de l'algorithme génétique en fonction du nombre de 
génération est illustrée ci-dessous : 
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gJ Genetic Algorithm 
File Edit View Insert Tools Desktop Window Help 


Best: 1.2165 Mean: 1.2179 


+ Best finess 
= Hean iness 


Fitness value 


ee O 40 50 
| Stop | | Pause | Generation 





- À la fin de processus d’itération la solution optimale qui représente le volume 
minimum du mur f (x) de l’équation 1 est: 





- Les méthodes empiriques et classiques pour la conception de mur de soutènement 
réalisent un volume de 4.71 m? tandis la méthode d'optimisation avec algorithme 
génétique réalise un volume de 2.85 m*. 

- Ce qui fait un écart de 1.86 M? en mètre linéaire, l'étude a été faite en Algérie et 
selon leur marché algérien, le mètre cube du béton armé coute jusqu’à 40000 DA. 

- Le tableau suivant montre le montant du gain qu’on peut gagner en utilisant AG 
dans l'optimisation de cet ouvrage si la longueur du mur est de (130 mètre linéaire) 








Volume du Prix unitaire | Montant (DA) | Montant 
mur (MÌ ) (DA) global 
Bureau Ai 40000 188400,00 24492000,00 
d'étude 
Algorithme 2.85 40000 114000 ,00 14820000,00 
génétique 
Gain pour 130 | 241.8 40000 9672000,00 9672000,00 


ml 


Tableau comparative 


Conclusion générale : 


En guise de conclusion l'optimisation des constructions en génie civil est un 
domaine d'actualité et de la recherche scientifique qui permet d'étudier des 
problèmes à l’aide des fonctions objectif et les résoudre en utilisant des méthodes 
algorithmiques pour réduire le coût de fabrication d’une structure. 

Pour notre cas la conception de mur de soutènement cantilever a cout optimal est 
un problème géotechnique et assez laborieux avec le calcul complet et le 
dimensionnement et toutes ses vérifications, tandis que l'application des méthodes 
d'optimisation assure un gain de temps et de cout en localisant rapidement la zone 
proche de la solution optimale. 

Pour arriver à la conception sécurisée et optimale, c'était indispensable de faire 
un couplage de la théorie géotechnique de Rankine avec l'algorithme génétique, ce 
dernier permet de converger rapidement vers la zone privilégiée de l’espace de 
recherche, en vérifiant les facteurs de sécurité de l'ouvrage. 

Ensuite L'étude comparative était assez intéressante et elle a bien éclaircir le but 
de notre étude en optimisant un volume de 1.8 m? et le gain financier qu’on a pu 
atteindre, donc l'AG a prouvé que c’est l’un des meilleurs méthodes d'intelligence 
artificielle efficace dans la conception et peut être une alternative prometteuse. 


Élargissement de sujet : 


Le temps de calcul est très important, et souvent dans les projets de génie civil on 
préfère adopter des durées moins longues dans la conception pour éviter le retard et 
les pénalités et l’AG examine seulement une petite partie de l’espace de recherche, et 
l'identification de l’optimum global requiert peu de temps, ainsi le réglage des 
paramètres peut être réalisable en théorie mais irréalisable en pratique,d’ouù l'idée 
d'hybridation de notre algorithme et le coupler par d’autres algorithmes 
d'optimisation par exemple : essaim de particules (PSO) ou algorithme de Piyavskii qui 
sont plus rapide par rapport à l'AG comme on peut améliorer notre fonction objectif 
d’où une minimisation de temps de calcul. 


Application au Maroc : 


Au Maroc le Ministère d'Équipement et de Transport et l'office national des 
chemins de fer sont chargés des projets routiers et ferroviaires, et absolument les 
routes régionales ou nationales contient des zones montagneuses, donc la nécessite 
des murs de soutènement est indispensable pour sécuriser le chemin routier, et 
souvent le budget annuel consacré au ministère ou à l'ONCF ne sont pas assez 
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important, et l'obligation de faire des minimisations dans le montant global des projets 
parait primordial, la conception des murs cantilever le long d’une route est toujours 
existante et parfois très abondant si le relief est montagneux. D’après notre étude, 
l'AG a pu optimiser un volume de béton très important dans un tronçon d’une route, 
et d’après le marché marocain un mètre cube du béton armé coute jusqu’à 250 DH 
donc on peut obtenir un gain financier facilement grâce à l'optimisation. 
Malheureusement au sein des bureaux d’études au Maroc existe rarement et sauf 
dans les organisations multinationales comme : NOVEC, CID. 


Donc l'application des méthode d'optimisation par tous les bureaux d’études sera un 


vif succès au niveau temps et cout et permet un dimensionnement qui constitue un 
problème de première importance pour l'ingénieur. 
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Annexe : 





La conception de murs de soutènement est basée sur le règlement parasismique (RPS 
2000), avant d'optimiser le volume de l'ouvrage la vérification des conditions liées à la 


stabilité, capacité et géométrie de la paroi est indispensable pour construire le mur et 
sont de nombre de trois : 


FSo > FSo design 
FSs > FSs design 
FSB > FSB design 


i) — FSo : Facteur de sécurité vis-à-vis du renversement : 

> Mr 

> Mo 

% Tel que : > Mr représente la somme des moments résistants et > Mo la somme des 
moments déstabilisants 


FSo = 


ii) — FSs : Facteur de sécurité vis-à-vis du glissement : 
> FR 
LI; 


% Tel que > Fp représente la somme des forces horizontales résistantes et > F) la somme des 
forces horizontales déstabilisantes : 


FSs = 





A 2.B.C 
X Fr = È Wwan). tan (ee) ++ P, ; E Fp = Pa- cosp 


% Où > Way est le poids total du mur, Øpase est l'angle de frottement du sol de fondation, B 
est la largeur totale de la semelle, C représente la cohésion de l'interface semelle et le sol de 
fondation, P, est la force de butée tandis que P, représente la force de poussée. 


iii) — FSB : Facteur de sécurité vis-à-vis de la portance du sol : 


FSB = Tu 
Imax 
% Où q, est la capacité portante ultime du sol de fondation et qmax est la pression 


maximale générée due aux forces agissantes sur le mur de soutènement et sont 


déterminée à l'aide des équations ci-dessous : 


YV 6.e 
dnin = UE) 


XV 6.e 
max 7 ra. + a) 
— ÿ V: La somme des forces verticales (résultante du poids du mur, le sol au-dessus de la 


Avec : 


semelle et la surcharge Q) 
—B: La largeur de la base ps B $} Mp- È Mo 


—e: L’excentricité du système de force résultant : 2 2 V 


Où, FSo design, FSs design,FSB design des facteurs de sécurité minimales qu’on doit 
absolument les dépasser pour assurer la stabilité de l'ouvrage. 
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